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Abstract 

We investigate the background dynamics when dark energy is coupled to dark matter in the 
universe described by Einstein cosmology and Loop Quantum Cosmology. We introduce a new 
general form of dark sector coupling, which presents us a more complicated dynamical phase space. 
Differences in the phase space in obtaining the accelerated scaling attractor in Einstein cosmology 
and Loop Quantum Cosmology are disclosed. 
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I. INTRODUCTION 

Our universe is undoubtedly undergoing an accelerated expansion driven by a yet unknown dark energy 
(DE)[l|, 0, 0, 4|. This mysterious energy component occupies almost 70% of the content of the universe today. 
The leading interpretation of such a DE is a cosmological constant with equation of state (EoS) u> = — 1. 
Although this interpretation is consistent with observational data, at the fundamental level it fails to be 
convincing. The vacuum energy density is far below the value predicted by any sensible quantum field theory, 
and it suffers the coincidence problem, namely, "why are the vacuum and matter energy densities of precisely 
the same order today?" . To overcome the coincidence problem, some sophisticated dynamical DE models 
relating the DE to scalar fields have been put forward to replace the cosmological constant [5[. 

Considering that DE contributes a significant fraction of the content of the universe, it is natural to look 
into its interaction with the remaining fields of the Standard Model in the framework of field theory. The 



possibility that DE and dark matter (DM) can interact has got growing attention recently [6] 



2l| . It has been 



argued that an appropriate interaction between DE and DM can influence the perturbation dynamics and 



affect the lowest multipoles of the CMB spectrum 9 



12} . Recently, it has been shown that such a coupling 



can be inferred from the expansion history of the Universe, as manifested in the supernova data together with 
CMB and large-scale structure [13(. Si gna tures of the interaction between DE and DM in the dynamics of 



galaxy clusters has also been analyzed 



14 Il5j | . It has been argued that the coupling between DE and DM 



can provide a mechanism to alleviate the coincidence problem and lead to an accelerated scaling attractor 



solution with similar energy densities in the dark sectors today 7 



lit 



A general interaction between DE and DM can be described in the background by the balance equations 



p x + 3H(1 + Lo x )p x = -Y, 

p m + 3Hp m = T, 



(1) 



where p m and p x correspond to the energy densities of DM and DE, respectively. Here T describes the coupling 
between DE and DM. Since the nature of dark sectors remain unknown, there is as yet no basis in fundamental 
theory for a specific coupling in the dark sectors. All coupling models discussed at the present moment are 
necessarily phenomenological 6|. There are two criterions to determine whether some models can be more 
physical justification than the others. One is to confront observations. The other is to examine whether the 
coupling can lead to accelerated scaling attractor solutions 16( , which is a decisive way to achieve similar energy 



densities in dark sectors and alleviate the coincidence problem. In this work we introduce a new form of dark 
sector coupling, T — 3cH p" pl^ a . This model is more general than the coupling discussed in the literatures. 
When a = 0, 1, it reduces to the cases with coupling between dark sectors solely proportional to the energy 
densities of DM 

12011 - 



18|, respectively. These two limiting cases have been examined thorough! 

3 



against observations and their possibilities to alleviate coincidence problem have also been discussed 
Here we will investigate the background dynamics when the DE is coupled to DM via this general interaction 
form. We will show that the general coupling leads to a more complicated dynamical phase space. 

Besides the discussion of the dynamics of DE with our general coupling to DM in the universe described 
by the Einstein theory, we will also extend our investigation to the Loop Quantum Cosmology (LQC). The 



LQC 



22 



23 



25 



2G 



2jJ in the cosmological context, 



24j | is the application of the Loop Quantum Gravity 
which keeps the properties of non-perturbative and background independent quantization of gravity. Recent 
investigations have shown that the loop quantum effects can be very well described by an effective modified 
Friedmann dynamics. There are two types of modification to the Friedmann equation. The first one is based 
on the modification to the behavior of inverse scale factor operator below a critical scale factor a» . Considering 
these modifications one can obtain many interesting results including the replacement of the classical big ban 
by a quantum bounce with desirable features 



29j . easier inflation 



ad, 



avoidance of many singularities 
and so on. However, the first type of modification to Friedmann equation suffers from gauge dependence which 
can not be cured and thus yields unphysical effects. The second type of modification to Friedmann equation 
is discovered very recently. It adds a —p 2 /p c term in the standard Friedmann equation which essentially 



24 



31 



32| . When energy density of the universe 



encodes the discrete quantum geometric nature of spacetime 
becomes of the same order of a critical density p c , this modification becomes dominant and the universe begins 
to bounce and then oscillates forever. Thus the bi g ba ng singularit y, t he big rip and other future singularities 



at semi-classical regime can be avoided in LQC 



24 



31 
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35 



. Therefore by using the second type 



of modification to Friedmann equation, the physically appealing features of the first type are retained. For 
the universe with a large scale factor, the first type of modification to the effective Friedmann equation can 
be neglected and only the second type of modification is important. Thus the dynamics of DE (phantom 



35 



37 



38, 39], quintom and hessence [40]) have been investigated recently in LQC on the basis of the second 
type of modification. It finds that the dynamical properties of dark energy models in LQC have different 
behaviors from those in the classical Einstein cosmology. Here we will examine the background dynamics of 



the LQC dominated by DE and DM where there is the general coupling L between dark sectors and compare 
the results with those in Einstein cosmology. 

The paper is organized as follows: in sections II and III, we study the dynamics of the interacting dark 
energy model in Einstein Cosmology and the LQC, respectively. In Sec. IV, we present numerical pictures 
of dynamics in Einstein cosmology and LQC. Our conclusions and discussions will be presented in the last 
section. 



II. DYNAMICS OF THE INTERACTING DARK ENERGY MODEL IN EINSTEIN 

COSMOLOGY 



In the Einstein theory, the universe is described by the standard Fricdmann equation 



H 2 = -rp, 



(2) 



where H is the Hubble parameter, p = p rn + p x is the total energy density and the constant k = 8ttG. Since 
we are concentrating on the late time accelerating universe, we have neglected the radiation and baryons for 
simplicity. 

Differentiating Eq.@ and using the conservation law of the total energy p + 'iH(p +p) = 0, we have 



H 



;{P + P)- 



(3) 



To analyze the evolution of the dynamical system, we introduce the dimensionless variables 



/np x 



d 



1 d 



(4) 



' V3H ' dN HdV 

where N = In a is the number of e-folding to represent the cosmological time. Using the above definitions, the 
Hubble equations can be rewritten as 



u 2 + v 2 = 1, 



(5) 



and 



H 

The effective total EOS u>tot is given by 



UJ X U 



U + V 



--(l+u x u 2 ). 



(G) 



Utot 



U x p x 



Px + Pm U 2 + V 



U) X U 2 
UJ X U . 



(7) 



Using dimensionless variables, the dynamical equations of the system can be expressed as 



3k 

T 

3v 

T 



L0 X (U 2 - 1) - C 



2 \ a— l-i 



LU X U + C[ 

V 



2 \ a 



(8) 



where the prime denotes a derivative with respect to N. The critical points u c , v c satisfy v! — and v' = 0. In 
order to study the stability of the critical points, we expand about the critical points u = u c + Su, v — v c + Sv 
and linearize the above equations near the critical points so that 



Su' 
Sv' 



3uj x ii 2 c - 1- (2a- l)c — 



2a-2 



Su 



3c(a-l)[ — 



2a-l 



Sv, 



3uj x u c v c + 3ca 



Su 



io x u\ — (2a — l)c 



5v. 



(9) 



The eigenvalues of the matrix of coefficients of the above equations encode the behavior of the dynamical 
system near the critical points. 

In general, for an arbitrary a it is difficult to obtain the analytical forms of the critical points. Here we 
only consider some specific values of a and examine the dynamics of DE with an interaction with DM. 

A. Case I: a — 



When a = 0, the interaction form reduces to the simple coupling between DE and DM in proportional to 
the energy density of DM. This simple interaction form has beenconfronted to observations and its possibility 
to alleviate the coincidence problem has also been examined 
can obtain two critical points for the specific model, namely: 



19(. From the dynamical system equations, we 



• Point A : (1, 0) 

• Point Bq : 



c I c 

LU X V Wsi 



(10) 



The eigenvalues of the coefficient matrix of the linearized equations around these critical points can be ex- 
pressed respectively as 



• Point Aq : Ai = 3w x , A2 = -(c + u> x ), 



• Point Bn 



Ax = -3c, A 2 = — 3(c + oj x ). 



(11) 



For point Aq, when c < — u> x , both eigenvalues Ai and A2 are negative, which indicates that Aq is a stable 
point. From Eq.([7]), we learn that the effective total EOS at point Aq is uitot = u x . Therefore we obtain 



6 

a cx — (1 + 3u> x ) £3(T+^7~ 2 a nd p cx a~ 3 ( 1+Wj; ). It means that point Aq is an accelerated scaling solution as 
u) x < —1/3 and there is singularity in the finite future as u> x < — 1. When c > the sign of Ai is always 

opposite to the sign of A2, which leads Aq to a saddle point. For point Bq, the critical point v c exists only 
provided that c < —uj x , which leads the sign of Ai always opposite to that of A2. Thus point Bq is a saddle 
point. 

B. Case II: a = 1 

In this limiting case the coupling between dark sectors is in proportional to the energy density of DE. This 
simple coupling has been examined using observational data and its effect to alleviate the coincidence problem 
has been discussed [191 ] . One can obtain two critical points in this dynamical system: 

• Point At : (0, 1), 



• PointB!: ( J1 + -, J-— ). (12) 



The eigenvalues of the coefficient matrix of the linearized equations are 

3 

• Point Ai : Ai = 0, A2 = — — (c + u) x ), 

• Point Br : Ai = 3(c + w B ), A 2 = 3(c + u; x ). (13) 

It is easy to examine that the critical point A\ is not a stable point, while B\ is stable when c < -Wj. The 
total effective EOS at point Bi reads uj to t = c + lu x . When lo x < —1/3 — c, point Bi can be an accelerated 

scaling solution. From Eq.© we find H — -rrr, rrn — jti H = —^rr* rrn — ttt- Thus as lo x < — 1 — c the 

universe will undergo super-accelerated expansion (H > 0) and end in the big rip. 

C. Case III: a = | 

Solving equations v! = and u' = 0, one can obtain two critical points of the dynamical system: 
• Point A 2 : 




Through analysis of the eigenvalues of the coefficient matrix, we find that point A2 is stable when c < — ^ , 



while point Bi is a saddle point. At point A2, the total effective EOS is uJtot — (w x — ^Jlli\ — 4c 2 )/2, which 



shows that if c > h and u> x < — 2c or c < | and — | < < — i — 3c 2 , we have LOtot < — \ so that A is an 



accelerated scaling attractor. When ui x < |(1 — 2^1 + 3c 2 ), one obtains w tot < —1, H > and finds there is 
a future singularity in this case. 



D. Case IV: a 



In this case the system has two critical points: 



• Point A 3 : (1, 0), 

• Point B 3 : 



' c 



2 ~ v 'c? — Au> x 



Auj t c 



(15) 



1io x y 2w x 

In term of the signs of the eigenvalues of the coefficient matrix, we find that point ^3 is stable when io x < 0, 
however, point B 3 is a saddle point. From Eq.([7]), we learn that point A3 can be an accelerated scaling 
attractor provided that uj x < — | and there is a future singularity asw x < — 1 . 

E. Case IV: a = 2 



When a = 2, the system has two critical points: 



• Point At : (0, 1) 

• Point B4 : 



1 - 



c — \fc 



-Vc 2 - 4uj x 



(16) 



2ui x V 2u) x 

Similarly, we find that point A 4 is unstable, while point B4 is stable for all uj x < 0. The total effective EOS 



at point B4 is u>tot = Wj — (c — \/c 2 — 4uj x c)/2. When ui x < — i(l + \/3c), we have w tot < — | and i?4 can be 
an accelerated scaling solution. Moreover, we find that u>tot > —1 as > —1 — and there is no future 
singularity in this case. But as lu x < — 1 — \fc one can obtain Lo to t < — 1 and a future singularity is inevitable. 
Thus in the Einstein cosmology the presence of the coupling terms can not remove the singularity entirely. 

III. DYNAMICS OF THE INTERACTING DARK ENERGY MODEL IN LQC 



In this section we are goin 



Friedmann equation into 



24, 



? t o extend our discussion to the LQC. The loop quantum effect modifies the 



31 



32 



33] 



3 V Pc 



(17) 



where p c = 167T ^ G2h is the critical loop quantum density and 7 is the Barbero-Immirzi parameter. Let us 
note here it has been suggested that 7 ~ 0.2375 by the black hole thermodynamics in LQG 34 1. 

Differentiating Eq. (fT"7j) and using the conservation equation of the total energy density p + 3H(p +p) = 0, 
where p = p x + p m , one can obtain 



H = --{p + p)[l- 2 -P 



(18) 



Adopting the dimensionless variables defined in (01, the evolution of the Hubble parameter in LQC becomes 



(u' + v z )[ 1 - 3H 



Pc 



= 1, 



and 



UJ X U 



U 2 + V 2 



The total effective EOS uJtot in LQC is given by 



(19) 



(20) 



Utot 



The dynamical system can be expressed as 



3« 
3v 

Y 



1 - U) x - c - 
\v 

la 



- 1 + C 



U) X U 
U 2 + V 2 



[(1 + lo x )u 2 + v 2 }{2 - u 2 - v 2 ) 

U 2 + V 2 



[(1 + uj x )u 2 + v 2 }(2 - u 2 - v 2 ) 



(21) 



(22) 



u 2 + v 2 

Obtaining the critical points and linearizing the system near them, we can study the stability of critical 
points by analyzing the first-order differential equations 



Su' = - 
2 



Sv' = 3 



2 - ' 2 ^ 2 + 3 ^-(2a-l) C p 

1 v 



— 3(uj x + l)u — V + 1 — LU X + 
2u> x u 3 v 



2Y2 



(u 2 +v 2 ) 2 



(a-l)c - 

v 



(u 2 + V 2 ) 
2a-l 

5v, 



2a-2n 



2uj x UV 3 

(lo x + l)uv + --5— — ^ + ac 

(u 2 + v i y 



2a-l 



I - - u 2 - 3v 2 + g „\ oxo y - (2a - l)c( - 



2<i 



(23) 



(w 2 + ?J 2 ) 2 

Solving the eigenvalues of the coefficient matrix of the above equations, we can know the behavior of the 
dynamical system near the critical points. Comparing with the Einstein theory, we find that the same critical 
points will have different eigenvalues in the coefficient matrix in LQC. This means that the dynamical property 
of the system in the LQC is different from that in the Einstein cosmology. 



9 

As did in section II, here we will also focus on some specific cases for simplicity, such as a — 0, 1, ^, — 1, 2. 
For these specific cases, the critical points in the LQC are the same as those in the Einstein cosmology. 

A. Case I: a — 

In this limiting case, the critical points of the dynamical system are 

• Point A Q : (1, 0), 



Ul x V U> x 



• Point B : ( J- — , Jl + — )• (24) 



The eigenvalues of the linearized equation around the critical points read 

3 

• Point A : Ai = -3(l + w x ), X 2 = -(c + lj x ), 

• Point B Q : Ai=3(c-1), A 2 = -3(c + ui x ). (25) 

It is easy to see that for the DE with ui x > — 1 point A is stable if c < -w x . The effective total EOS uu to t = w x , 
thus the critical point A is an accelerated scaling solution without a future singularity for — 1 < u x < — 1/3. 
When the DE is phantom like u) x < — 1, Aq is a saddle point. This is different from that in the Einstein theory, 
where A is still stable for the phantom DE provided that c < —uj x . Point B is a saddle point if c < 1 and 
is unstable when c > 1. In the Einstein theory, B is always a saddle points. This means that the presence 
of the term — in the Fricdmann equation due to the quantum correction changes the dynamical properties 
of autonomous system. Moreover, when the value of u x goes beyond the ranges of — 1 < w x < —1/3, the 
quantum bounce originated from the term — ^ will leads to the avoidance of the future singularity. 

B. Case II: a = 1 

The critical points are the same as those in the Einstein theory 

• Point A x : (0, 1), 



U! x V UJ 



• Point Bi : ( Jl + — , )• (26) 



However, the eigenvalues of the coefficient matrix of the linearized equations near the critical points become 

3 

• Point A x : Ai = -3, A 2 = --(c + w x ), 

• Point Bi : Ai = 3(c + u x ), A 2 = -3(1 + c + u x ). (27) 



10 

Point Ai is a stable point when c > — u> x . From Eq. (j20[) we have io to t — at the point A±, which means 
that in this model our universe described by LQC will enter DM dominated era and there is no singularity 
in the finite future. However, in the Einstein theory A\ is unstable. Point B\ can be stable provided that 
— (1 + u> x ) < c < —lj x . This region of c to keep B\ stable is smaller than that in the Einstein theory. At point 
Bi, we have ui to t = c + u> x> which shows that when — 1 — c < u> x < — -| — c, we have — 1 < oj to t < — h, so that 
point B\ corresponds to an accelerated attractor without a future singularity. 

C. Case III: a = \ 



The critical points can be found at 



Point A 2 



• Point B 2 



\ 



11/ 4c 2 



1_1 



\ 



11/ 4c 2 



1 1 



4c 2 



(28) 



Analyzing the stability we find that when DE is of quintessence type, point A 2 is stable provided that c < 



When the DE is of phantom type, A 2 can be stable only when + oj x ) < c < — and ui x > —2. From the 

total effective EOS, we learn that when c > | and — 1 — c 2 < u x < — 2c or c < | and — 1 — c 2 < u> x < — | — 3c 2 , 
point ^2 is an accelerated attractor. Similarly, when the value of w x goes beyond the ranges above, the effects 
from the term —y- will make the future singularity disappear. Point B 2 is a saddle point, which agrees with 
that found in the Einstein theory. 



D. Case IV: a = -1 



Critical points of the system read 



• Point A 3 : (1, 0), 

• Point B 3 : 



' c 



1 



(29) 



2oj x y 2uj x 

For the quintessence type DE, point A 3 is a stable point. Since uitot — Wj, ^3 is an accelerated scaling 



attractor. For the phantom type DE, A 3 is a saddle point when c < — t-^ — and is unstable when c > — . 1 
Point i?3 is a saddle point always. 
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E. Case IV: a = 2 



When a = 2, the system has two critical points: 
• Point A 4 : (0, 1), 



Similarly, we find that point A is a saddle point as that in the Einstein theory. When the DE is of quintessence 
type, point B4 is stable and when — I < lu x < — 1(1 + V3c) B4 is an accelerated scaling attractor. When the 
DE is of phantom type, B4 can be stable only when c > (1 + uj x ) 2 . When —(1 + y/c) < ui x < — 1(1 + \/3c), 
B 4 can be an accelerated scaling solution as well. Although oj to t < -1 aswj < -1- \fc, the loop quantum 
effects will cancel off the future singularity in the evolution of the universe. Thus in the LQC the big rip can 
be removed entirely. 

IV. NUMERICAL RESULTS 

In this section we confirm numerically the complicated stability conditions for critical points obtained above. 
For the new general form of the interaction between DE and DM, T = 'icH p x p]^ a > we have more complicated 
dynamical phase spaces. We find that the position of the critical point and its stability depend not only on the 
coupling constant c, the DE EOS lo x and the exponent a in the coupling, but also on the theory to describe 
the universe. In figure (1), we show the stable regions in the parameter space (c, u x ) by choosing a = 0.5 and 
2. In Einstein cosmology the critical points A2 (a = 0.5) and B4 (a = 2) are late time attractor in the region 
I + II. However in LQC, A 2 and B4 are late time attractors only in the region II. In LQC we see that the 
region of the location of the accelerated scaling attractor has been reduced compared to the Einstein theory. 
This holds true for other values of a. 

In figure (2), we plot the numerical result to illustrate the phase space trajectories for our coupling model 
with chosen c and uj x in the stable region and selected a in the Einstein cosmology and the LQC. Since the 
LQC has the same critical points as that in the Einstein cosmology, we see that for fixed a, all trajectories with 
different initial conditions in the Einstein cosmology and LQC converge to the same final state determined by 
parameters c and uj x . This means that our universe will enter an era with similar energy densities of DE and 
DM. 

In figure (3), we show the evolution of the effective EOS u> to t — u x p x /(p m + p x ) in the stable region. We 



12 

see that in the final state the total state parameter uj tot tends to a constant, which is determined by values 
of 

Px c i Pm c , u x and a. For selected values of u) x ,c to be within a certain range discussed in sections II and 
III for different values of a, such as u x = —0.6 and c = 0.18, we can have stable critical point and meanwhile 
we can get the total effective EOS LOtot < -| in the final state as displayed in figure (3). Our universe will 
enter a final state with a constant energy ratio between DE and DM and accelerate forever for all chosen as. 
However when values of to Xl c are beyond the range discussed in Einstein cosmology and LQC, our universe 
will enter a decelerated expansion. 

In figure (4) we exhibit the evolution of p(t) for chosen a = — 1, p c = 0.82 and parameters (uj x ,c) in 
the unstable region in LQC. We find that the universe finally enters an oscillating regime in the LQC. The 
oscillating frequencies of p(t) depend on the coupling constant c and DE EOS lo x . This oscillation behavior 
makes the universe experience bouncing, which can avoid singularity faced in the usual Einstein cosmology. 
This property also holds for other values of a in LQC. 

V. CONCLUSIONS AND DISCUSSIONS 

In this paper we have studied the background dynamics when DE is modelled coupling with DM via a new 
general form 2>cH p x p]^ a in Einstein cosmology and LQC. For selected values of a, we have examined stability 
behaviors of critical points and found accelerated scaling solutions to account for the similar energy densities 
in dark sectors today. In LQC, the parameter space for the existence of the accelerated scaling attractor is 
found smaller than that in the Einstein cosmology. In the unstable region, the universe described by the LQC 
will enter an oscillatory regime which can help to avoid the singularity usually met in Einstein cosmology. 

The background dynamics for the new general form of dark sector coupling leads to more complicated 
features in dynamical phase space. In order to confront this model with observations, the cosmological 
perturbations with this coupling form need to be disclosed. Some discussions in this direction has been 
addressed recently in 



41|. 
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FIG. 1: The stable regions in the (ui x , c) parameter space for fixed a (the left for a = 0.5 and the right for a — 2). In 
Einstein cosmology the critical points A2 (a = 0.5) and B4 (a = 2) are late time attractor in the region I+II. But in 
LQC, A2 (a = 0.5) and B4 are late time attractor only in the region II. The region III represents the region of the 
solution without physical meaning. 
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FIG. 3: The evolution of total cosmic energy to for uj x = —0.6, c = 0.18 and fixed a (the left for Einstein cosmology, 
and the right for LQC). 
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